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Abstract:  We  use  a  simple  model  to  study  the  color  change  taking  place  when 
sodium  atoms  are  absorbed  in  the  zeolite  sodalite.  The  Hamiltonian  is  that  of 
an  electron  moving  in  the  electrostatic  field  created  by  the  ions  in  the  zeolite 
firamework  and  by  the  alkali  ion  core.  We  examine  to  what  extent  the 
absorption  spectrum  is  sensitive  to  the  magnitude  of  firamework  charges,  the 
orientation  of  the  Na^  cluster  in  the  sodalite  cells,  the  localization  of  the 
electron,  the  nature  of  the  alkali  impurity  (Li,  Na,  K),  and  the  laser 
polarization.  Comparing  the  calculated  spectra  with  the  experimental  ones 
helps  decide  which  framework  charge  models  are  consistent  with  the 
absorption  spectrum. 
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L  Introduction 

Sodalite  is  a  zeolite  with  the  stoichiometric  formula  (Na3)^'^  ([AlSi04]*)3. 
It  consists  of  a  periodic  arrangement  of  cages  having  the  structure  shown  in 
Fig.  1.  When  sodalite  is  exposed  to  sodium  vapor  some  of  the  cages  trap  a 
soditim  atom  to  form  the  compound  (Na^l'^e*  ([AlSi04]‘)3-  Here  e*  denotes  an 
electron  that  is  no  longer  tightly  boxind  to  the  absorbed  sodium  atom.  In  what 
follows  we  assume  that  the  (Na4)^'*‘e*  cluster  in  this  compound  is  a  tetrahedron 
whose  center  coinddes  with  that  of  an  alumino-silicate  cage.  This  assumption 
is  based  on  the  analogy  with  the  compound  hydroxysodalite,  (Na4)^^(OH)* 
[AlSi04*]3,  which  has  the  same  formula  as  the  Na  doped  sodalite,  with  an  OH’ 
group  instead  of  an  electron.  X-ray  diffraction  studies^^^  of  (Na4)^'^(OH)’ 
CAlSi04*33  show  that  the  cx3rgen  atom  of  the  (Na4)^OH*  duster  is  in  the  middle 
of  the  sodalite  cage,  and  the  Na  atoms  are  at  the  comers  of  a  tetrahedron 
surrounding  the  OH*  group.  ESR  studies®^  of  (Na4)^‘'’e’([AlSi04]*)3  show  that  e’ 
is  located  at  the  center  of  the  (Na4)^'‘'  duster,  just  like  OH’  in  (Na4)^'^(OH)’. 
These  similarities  prompt  us  to  asstime  that  the  ions  in  (Na4}^‘*'e’  have  the 
same  tetrahedral  structure  as  those  of  (Na^)**  (OH)*. 

Doping  sodalite  with  sodium  causes  a  change  of  color,  from  transparent 
when  sodium  is  absent,  to  light  blue  when  the  doping  level  is  low,  and  to  black 
when  the  sodalite  is  saturated^^’^^  The  process  is  reminiscent  of  the  ionization 
of  sodium  dissolved  in  liquid  ammonia,  which  turns  blue  because  the  solvated 
electron  absorbs  light. 
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Recently  Srdanov  et  have  measured  the  absorption  spectrum  of  Na 
doped  sodalite  in  the  low  concentration  limit.  In  this  article  we  use  a  simple 
model  for  the  electron  -  sodalite  interactions  to  calculate  this  spectrum.  The 
calculations  have  several  goals.  First,  we  want  to  understand  the  spectrum  in 
terms  of  the  electron  dynamics  in  zeolite.  Second,  we  wish  to  find  which  sets  of 
framework  charges  are  consistent  with  the  observed  spectrum.  Third,  we 
would  like  to  know  how  sensitive  the  spectrum  is  to  the  structure  of  the  Na4 
cluster.  Fourth,  we  wish  to  know  whether  the  polarization  dependence  of  the 
absorption  cross  section  contains  enough  information  to  make  absorption 
experiments  with  single  crystal  zeolite  worthwhile. 

The  properties  of  an  electron  "dissolved**  in  a  cltister  are  of  great  current 
interestf^^.  The  theoretical  work  on  this  subjectf*^  is  rather  sophisticated,  but 
uses  several  approximations  to  make  the  calculations  feasible:  the  potentials 
are  schematic,  the  adiabatic  approximation  is  made  and  the  nudear  degrees 
of  freedom  are  sometimes  treated  classically.  To  find  whether  the  results  of 
such  simulations  are  satisfactory  it  is  necessary  to  compare  them  to 
experiments.  Unfortunately,  it  seems  to  be  difficult  to  achieve,  in  gas  phase 
systems,  the  cluster  density  needed  for  the  spectroscopic  measurements. 
Furthermore,  one  is  not  qmte  certain  that  the  clusters  prepared 
experimentally  have  the  same  structure  as  the  ones  studied  theoretically. 

In  view  of  these  difficulties,  the  clusters  encountered  in  zeolite  systems 
have  some  advantages:  it  is  not  difficult  to  prepare  samples  with  high  cluster 
densities;  moreover,  since  the  clusters  are  likely  to  have  a  well  defined 
orientation  the  laser  polarization  can  be  used  to  increase  the  dynamic 
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information.  In  gas  phase  experiments  this  information  is  degraded  by 
orientational  averaging.  Finally  X-ray  difEraction  measurements  can 
determine  the  structure  of  the  cluster  and  its  orientation  with  respect  to  the 
framework. 

The  doped  sodalite  discussed  here  is  a  material  of  interest  on  its  own. 
One  can  think  of  it  as  a  collection  of  (Na4)^e‘  clusters  embedded  in  a  quasi- 
rigid  ionic  "medium"  provided  by  the  alumino-silicate  framework.  The 
clusters  occupy  the  sites  of  a  periodic  lattice  and  their  concentration  can  be 
varied.  At  higher  concenfrations,  the  probability  that  two  clusters  occupy 
neighboring  lattice  sites  is  high.  The  electron  wave  functions  of  such 
neighbors  overlap  and  this  affects  the  optical  and  the  electron  transport 
properties  of  the  material.  The  geometric  and  electrostatic  parameters  in 
these  structures  are  very  different  from  those  encountered  in  metals  or 
semiconductors.  Moreover,  they  can  be  changed  by  modifying  the  shape  and 
the  size  of  the  cages,  and  the  charges  on  the  framework.  It  is  reasonable  to 
hope  that  such  materials  will  have  imusual  physical  properties. 

The  article  is  organized  as  follows.  In  Section  11  we  explain  the  theory 
used  to  calculate  the  absorption  spectrum  and  in  Section  HI  we  present  the 
model  used  for  the  electron  -  zeolite  interactions.  Some  computational  details 
are  discussed  in  Section  IV  and  the  restilts  of  the  calculations  in  Section  V. 

n.  THEORY 

We  consider  a  rigid  sodalite  framework  containing  an  excess  Na  atom. 


The  Hamiltonian 
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H*p^/(2m)  +  V(r;R,Q)  01.1) 

contains  the  kinetic  energy  of  the  electron  plus  its  interaction  energy  with  the 
framework  atoms.  The  electron  position  and  momentum  are  r  and  p, 
respectively.  The  positions  of  the  Na  ions  and  of  the  framework  atoms  are 
symbolically  indicated  by  R  and  their  charges  by  Q.. 

BecatJise  the  electron  excitation  gap  (about  2  eV)  is  much  latter  than  kgT  ( 
0.026  eV  at  300K),  the  electron  is  initially  in  the  ground  state  (x,y,2 1  g)sg(x,y,z). 
We  calculate  this  from 

<x,y,z  I  g>  a  lim  (x,y,2 1  exp  [-P  H]  I  'P)  (II.2) 

P-»oo 


where  p  is  a  real  number  with  units  of  energy*^.  For  (x,y,2 1 T)  s  'P(x,y,z)  one  can 
use  any  normalizable  function  that  ba^  a  nonzero  projection  on  the  groimd  state. 
The  numerical  procedure  used  to  evaluate  I  g)  from  (II.2)  has  been  described  by 
HeUsing,  Nitzan  and  Metiu^^^. 

To  calailate  the  absorption  cross  section  aioai)  we  use  a  slight  modification 
of  Heller's  formula^®^ 


dt  expCicot)  C(t) . 


(n.3) 


The  overlap  integral  C(t)  is: 
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at)*  exp(iEjt/h)<<D(0)  I 


(n.4) 


with 


1^,  0>*e*rlg> 


(n.5) 


and 


I  <&,  t  >  s  exp(-iH  t  ^  )  1 0, 0  >  . 


(n.6) 


Eg  is  the  ground  state  energy  of  the  electron  and  e  is  the  polarization  vector  of  the 
laser.  The  propagator  in  (II.6)  is  evaltiated  with  the  Fleck,  Morris  and  Feit^^^^ 
method.  We  call  1 0, 0 )  the  promoted  state. 

Heller  has  pointed  out^^^  the  connection  between  the  resolution  of  the 
absorption  spectrum  and  the  length  of  time  in  which  the  promoted  state  explores 
the  potential  energy  surface.  This  was  used  by  Engel,  Schinke,  Henning  and 
Metiu^^^^  to  obtain  dynamic  information  by  generating  a  sequence  of  low 
resolution  spectra: 


a(<o;  x)  -  CO  Re 


f 

k 


dt  £^(t)  exp(itot)  C(t). 


(n.7) 


The  "window"  function 
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£5(t)*exp[-(t/t)2].  (n.8) 

cuts  off  the  time  evolution  of  the  overlap  integral  at  times  longer  than  x.  In  the 
frequency  domain,  this  broadens  the  spectrum  and  gexmrates  a  version  having 
the  resolution  A(0s2/t.  Thus,  a  spectrum  taken  with  the  resolution  do  contams 
information  about  the  motion  of  the  promoted  state  finr  a  time  t  s  2/Aod. 

The  measurements^^^  are  made  on  a  sodalite  powder.  To  compare  our 
calculations  with  the  results  of  the  experiment  we  evaltiate 

N 

o(fl);  x)  s  (1/N)  ^  0^(0);  x ).  (II.9) 

nal 


The  cross  sections  corresponds  to  a  particular  iwlarization  direction  with 
respect  to  the  zeolite  fr'amework.  The  sum  is  over  a  representative  sample  of 
directions  and  Na3  orientations  (see  Section  V.6) . 

m.  THE  MODEL  FOR  V(  r ;  R,Q) 

The  electron  interaction  with  the  framework  and  the  alkali  ion  is 
described  by  two  body  spherically  symmetric,  energy  independent,  local 
pseudopotentials^^^. 

The  energy  of  the  interaction  between  the  electron  and  the  framework  ions 
is 


V(r;  {R,Q}  )=  ^  |  |  /xf  ]  /  minC  I  r-R^  ( ,  ].  (III.  1) 

i 
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The  atoms  forming  the  zeolite  frame  are  treated  as  point  charges 
located  at  the  positions  R^.  The  electron  charge  is  denoted  e  and  is  a  cutoff 
distance  which  truncates  the  coulomb  potential.  This  form,  suggested  by 
Shaw^^^^,  has  been  used  for  modeling  the  dynamics  and  absorption  spectrum  of 
an  excess  electron  in  a  variety  of  atomic  clusters^^^  The  cut  off  is  used  to 
eliminate  the  effect  of  remote  charges.  The  smooth  exponential  cut  off,  exp[  -( t  r- 
R^l/X)”  ],  is  preferable  when  the  numerical  method  uses  Fourier  transforms. 
We  use  X  s  15.0  -  25.0  A  and  n=4.  Including  the  zeolite  atoms  located  between  40 
and  53  A  from  the  grid  center  does  not  affect  the  spectrum.  To  calculate  the 
energy  of  the  electron«framework  interaction  we  use  about  16,000  zeolite  atoms; 
we  include  all  atoms  within  40A  from  the  grid  center. 

The  charges  Qj  and  the  cutoff  distances  R^^  used  in  Eq.dQ.!)  are  given  in 
Table  L  The  values  of  R^|  are  chosen  such  that  the  two  body  interaction  between 
the  electron  and  atom  "i"  approximates  the  ionization  potential  of  the 
corresponding  ion  (or  atom)  with  charge  Q|.  For  example,  if  we  deal  with  the  ion 
A^''’e*  the  cut-off  distance  is  chosen  so  that  the  interaction  between  e*  and  A^'*’ 
gives  a  reasonable  value  for  the  ionization  potential  of  A'*'.  The  values  chosen  for 
Na'*'  are  discussed  in  Section  V.l. 

Approximation  the  interaction  potentials  by  electrostatic  point  charge 
interactions  or  by  point  dipoles  goes  back  to  the  crystal  field  theory  (CFT)  of  Bethe 
in  the  1920's^^^^.  The  magnitude  of  the  charges  is  still  xmsettled.  Assuming  that 
the  ionic  charges  are  those  suggested  by  the  valence  of  the  atoms  (e.g.  3  for  Al,-2 
for  O)  fails  dramatically  for  some  cases  (as  it  would  for  ours  -  see  Section  V). 
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Usmg  his  definition  of  electronegativity,  Pauling^^^^  concluded  that  zeolites  are 
ionic  crystals  with  50%  ionidty.  Density  functional  theory  proposes^^^  a  diJSerent 
defizution  for  electronegativity  which  leads  to  different  charges.  Since  there  is  no 
widely  accepted  set  of  charges  we  calculate  the  absorption  spectrum  for  several 
charge  distributions,  denoted  by  a,  b,  c,  d,  and  e  in  Table  L  The  comparison  of  the 
computed  spectra  with  the  measured  ones  helps  establish  which  proposals  are 
compatible  with  the  observed  optical  properties. 

The  a-model  uses  the  smallest  charge^^®^.  These  charges  were  arrived  at 
by  regarding  the  sodalite  framework  as  a  neutral  Si02  firamework  in  which  some 
O  atoms  were  replaced  with  Al.  The  fiactional  charges  used  by  the  model  are 
caused  by  the  imbalance  generated  by  this  substitution.  The  alkali  counterions 
take  on  a  -fl  charge  to  maintain  charge  balance  with  the  fractional  negative 
charge  on  the  oxjrgens.  The  b-model  is  an  interpolation  between  the  a-model  and 
the  c*model.  The  c*model  is  based  on  the  work  of  Leherte,  et.  who  use 

MuUiken  ab  initio  STO-3G  atomic  net  charges  to  suggest  for  the  ferrierite  zeolite 
about  a  +1.5  charge  on  the  Si  atoms  and  a  *0.7  charge  on  oxygen.  The  d-model  is 
from  Van  Genechten,  et.  al.^^®^  who  use  electronegativity  equalization  methods  to 
place  a  +1.9  charge  on  the  Si  atoms  in  sodalite  and  a  -0.95  charge  on  oxygen.  The 
e-model  is  an  interpolation  between  the  d-model  and  a  model  suggested  by 
Skorczyk^^®^  who  places  a  +3.03  charge  on  the  Si  atoms  ,  +2.45  on  the  Al ,  and  -1.62 
on  the  oxygen. 

IV.  COMPUTATIONAL  DETAILS 

The  Feit,  Morris,  and  Fleck^^®^  split  operator  method  used  to  evaluate  how 
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exp(-iHt/h)  and  ezp(-PH)  modify  a  given  wave  function  use  a  time  step  of  0.0 Ifs 
and  a  ^  step  size  of  0.05  eV*^  .  The  three  dimensional  spatial  grid  has  48  points 
with  a  spacing  of  0.546  A  in  each  direction.  We  tested  the  convergence  of  the 
numerical  method  by  performing  calculations  with  smaller  and  smaller  time, 
and  spatial  steps  until  the  spectrum  did  not  change.  A  grid  with  64  points  and  a 
step  of  0.4095  A,  which  covers  3  unit  cells  in  each  4-fold  symmetry  direction,  and 
also  a  grid  covering  4  unit  cells  in  each  4-fold  symmetry  direction  shifted  the 
ground  state  energy  of  the  electron  by  0.02eV  (about  0.5%)  and  the  promoted  state 
energy  by  0.0  leV  (about  0.4%).  We  typically  propagate  the  system  for  a  time  length 
of  30  to  120  fs.  For  this  period  the  converged  results  conserve  energy  to  within 
about  0.2%. 

V.  RESULTS 

V.1  The  spectrum  of  the  isolated  Na  atom 

The  absorption  spectrum  of  the  isolated  Na  atom  is  obtained  by 
suppressing  in  the  Hamiltonian  all  but  the  Na  ion  -  electron  interaction.  We  have 
computed  it  to  test  whether  using  a  truncated  coulomb  potential  is  reasonable  and 
to  compare  it  with  the  spectrum  of  the  Na  atom  absorbed  in  sodalite. 

In  Fig.  2  we  show  the  absorption  spectrum  of  the  sodium  atom  with  the 
cut  off  parameter  R^|sl.725  A.  This  value  gives  an  ionization  potential  of  5.139  eV. 
The  time  constant  x  of  the  window  function  is  40  fs.  This  gives  a  line  width  on  the 
order  of  h/t. 

The  calcxilated  ionization  potential,  oscillator  strengths  and  peak 
frequencies  are  compared  to  the  measured^^^^  ones  in  Table  H.  The  oscillator 


strengths  and  frequencies  of  the  first  two  peaks  in  the  calculated  spectrum  agree 
with  the  measured  ones  quite  well. 

Since  the  pseudopotential  is  spherically  symmetric  the  dipole  selection 
rules  allow  transitions  frum  the  s  ground  state  to  the  np,  states  ( the  z  aria  was 
taken  parallel  to  the  laser  polarization).  The  transition  intensity  is  highest  for 
ls~»2p2  and  decreases  with  n. 

y2  The  dependence  of  the  absorption  spectrum  on  the  framework  charges. 

It  is  useful  to  think  of  the  sodium  doped  sodalite  (Na4)^''‘e*([AlSi04]*)3  as 
an  electron  in  an  ionic  system  consisting  of  a  tetrahedral  soditim  cluster  (Na^)^'*' 
imbedded  in  the  alumino  -  silicate  framework  of  the  zeolite.  To  understand  how 
the  framework  affects  the  electron  we  examine  first  the  isolated  (Na4)^''‘e'  system. 
The  dimensions  of  the  (Na4)‘^  tetrahedron  are  the  same  as  inside  the  zeolite.  This 
spectrum  serves  as  a  reference  used  to  examine  the  influence  of  the  zeolite 
framework.  The  absorption  spectrum  (Fig.  3)  of  was  calculated  for  the 

laser  polarization  shown  in  Fig.  4.  The  tetrahedral  field  of  the  ions  breaks  the 
spherical  symmetry  and  splits  the  peaks  appearing  in  the  spectrum  of  the  isolated 
Na  atom  (Fig.  2).  The  peak  positions  are  shifted  and  the  ionization  potential  is 
18.82  eV.  Considering  the  enormous  difference  in  the  ionic  charge,  the  similarity 
between  the  spectra  of  Na  and  that  of  e'(Na4)^''’  is  rather  surprising.  This 
similarity  is  partly  due  to  the  small  size  of  the  tetrahedron  which  makes  the 
potential  energy  felt  by  the  electron  almost  spherically  symmetric. 

Next  we  study  how  the  friamework  charges  affect  the  properties  of  the 
electron  when  the  (Na4)^''*e'  cluster  is  located  inside  the  sodalite.  The  spatial 
distribution  of  the  electron  in  the  ground  state  g(z,y,z)  is  characterized  by 
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Jf*’2 

r^drsin  (0)d6d<^  I  g(r,0,^)  I  (V.l) 

n 

The  values  of  and  T2  are  chosen  so  that  we  obtain  the  probability  of  finding  the 
electron  in  the  (Na4  tetrahedron  region,  on  the  nearest  alumino  -  silicate  shell 
or  in  the  region  outside  the  unit  cell  in  which  the  tetrahedron  is  located.  The 
results  of  these  calculations  are  given  in  Table  m. 

For  the  models  a  through  d  the  electron  tends  to  remain  localized  ixiside 
the  tetrahedron  (  see  Table  HE).  The  only  exception  is  model  e  for  which  the 
probability  that  the  electron  is  located  outside  the  unit  cell  containing  the 
tetrahedron  is  high.  ESH  ezperiments^^^  on  sodalite,  suggest  that  the  excess 
electron  is  localized  around  the  (Na4)^'^  tetrahedron.  Moreover,  the  absorption 
spectrum  for  the  model  e  is  rather  different  from  that  observed  experimentally. 
For  these  reasons  we  abandon  model  e  and  concentrate  in  what  follows  on  the 
models  (a-d). 

In  Fig.  5  we  show  the  spectra  of  the  doped  sodalite,  for  the  charge 
distributions  a,b,c,d  of  Table  I,  and  the  laser  polarization  shown  in  Fig.  3.  The 
magnitude  of  the  charges  is  progressively  higher  as  we  go  firom  model  a  to  model 
d  (Table  I).  The  spectrum  begins  shifting  dramatically  only  when  the  ionic 
character  of  the  zeolite  atoms  is  dose  to  50%  (model-d  of  Table  D.  Only  in  the  d- 
model  the  forces  due  to  the  framework  charges  are  large  enough  to  compete  with 
those  exerted  by  the  (Na4)'‘^  tetrahedron. 
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V3  The  absoiiption  spectrum  and  electron  dyxuanicstixne  scales 

We  examine  here  the  connection^^^^  between  the  absorption  spectrum  and 
the  time  scales  of  electron  dynamics. 

In  Fig.6  a,b,c  we  show  the  spectra  for  the  d-model  (see  Table  D  of  the 
doped  sodalite  for  three  different  resolution  tiznes,  x  s  5, 20,  and  40  is.  For  x  s  5  &, 
the  energy  resolution  is  on  the  order  of  1  eV  and  for  40  &  it  is  on  the  order  of  O.leV. 
The  presence  of  two  bands  in  the  spectrum  is  discerned  even  if  the  electron 
samples  the  potential  for  5  fs.  Sampling  for  20  fs  is  enough  to  resolve  the  presence 
of  three  peaks,  and  a  40  fs  sampling  resolves  all  the  structure.  Further  sampling 
will  only  narrow  the  lines  obtained  at  40  fs.  On  this  time  scale  the  nuclear  motion 
does  not  have  time  to  affect  the  spectrum.  If  the  motion  of  the  firamework  is 
included  it  will  contribute  low  frequency  phonon  side  bands  which  have  little 
information  regarding  the  electron  dynaznics. 

To  get  a  better  feeling  for  the  kind  of  motion  that  generates  the  spectrum 
we  ezazziine  the  evolution  of  the  mean  radial  electron  position  (r)^  for  the  time 
dependent  promoted  state.  In  the  e'(Na4)^'^  duster  (r)^  undergoes  localized 
osdllations  with  a  maximtim  amplitude  of  about  0.4  A.  When  the  duster  is 
imbedded  in  the  the  d-model  sodalite  the  electron  performs  more  extensive 
excursiozis  zmder  the  influence  of  the  framework  charges  with  a  maximum 
amplitude  of  about  1.7  A. 

hi  Fig.  7  we  plot  the  probability  P(r]^,r2,'t)  that  the  electron  in  the  promoted 
state  is  located  between  r^^  and  r2  (see  Eq.  (V.l)),  for  the  d-model.  Here  r  denotes 
the  radial  distance  measured  from  the  center  of  the  uziit  cell  (which  coincides 
with  the  center  of  the  Na^*^  tetrahedron)  Initially  the  promoted  state,  like  the 
electron  ground  state,  is  localized  around  the  tetrahedron;  83%  of  the  electron 
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density  is  within  3.4  A  of  the  center  (each  is  located  at  2.40lA  from  the 
center).  There  is  a  16%  probability  that  the  electron  is  in  the  region  of  the  closest 
zeolite  shell  ( between  3.4  and  5.4  A).  Only  1%  of  the  electron  density  extends  into 
the  surrounding  cages. 

In  the  first  five  fs  50  to  80%  the  electron  density  moves  into  the 
surrounding  cages  at  the  expense  of  the  population  around  the  tetrahedron.  This 
suggests  that  the  excited  states  that  cany  oscillator  strength  are  delocalized  and 
extend  into  the  neighboring  zeolite  cells.  For  the  isolated  Na  atom  or  the  bare 
(Na4)'^'^a-  cluster  the  probability  density  of  the  promoted  state  changes  very  little  in 
time.  Another  way  to  interpret  this  calculation  is  that  if  excited  with  a  5-function 
pulse  the  electron  will  move  out  of  the  Na^***  trap  in  five  femtosecond.  This 
suggests  that  if  the  trap  concentration  is  low  the  material  will  photoconduct. 

y.5.  The  dependence  of  the  spectrum  on  the  Himgn«inTig  of  the  Na^  cluster. 

So  far  we  have  assumed  that  the  (Na4)^''’  cluster  in  the  sodalite  cage  has 
the  same  dimension  and  orientation  as  the  cluster  in  hydroxysodalite.  We 

have  done  this  out  of  necessity,  since  no  X-ray  studies  of  the  doped  sodalite  are 
available.  The  absorption  spectrum  is  very  sensitive  to  the  distance  d]^  between  the 
Na  ion  and  the  center  of  the  cell  (which  is  the  same  as  the  center  of  the 
tetrahedron)  (Fig.8).  The  ground  state  energy  also  changes  from  -4.594eV  for 
d]^s2.60lA  to  -6.1388eV  for  d]^=2.30lA.  This  sensitivity  to  the  size  of  the  tetrahedron 
makes  the  X-ray  study  of  sodium  doped  sodalite  essential  for  a  detailed 
understanding  of  the  absorption  spectrum.  It  also  suggests  that  materials 
differing  through  the  composition  of  the  tetrahedron  are  likely  to  have 
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substantially  different  properties. 

V£.  Hie  dependence  of  the  spectrum  on  the  orifintation  of  the  Na3cfaiste3cs. 

The  spectrum  may  al^  be  influenced  by  the  orientation,  the  size  and  the 
positions  of  the  (Na3)^'^  clusters  located  in  the  cells  surrounding  the  one  occupied 
by  the  (Na4)'^'*'  cluster.  Since  the  (Na3)^'*'  clusters  do  not  show  up  in  diffraction 
studies,  we  may  assume  that  their  orientation  is  random.  The  observed  spectrum 
is  then  the  average  over  all  these  orientations.  As  shown  by  Fig.  9  the  spectrum 
changes  significantly  with  the  orientation  of  of  the  (Na3)^''‘  clusters.  Since  the 
groimd  state  of  the  electron  is  very  well  localized  around  the  Na4  tetrahedron,  the 
dependence  of  the  spectrum  on  the  orientation  of  the  Na3  clusters  in  the  adjacent 
cells  is  another  indication  that  the  ezdted  states  are  delocalized. 

The  spectrum  is  also  sensitive  to  the  size  of  the  (Na3)^''‘  clusters  (Fig.  10) 
from  the  cells  surrounding  the  one  in  which  the  tetrahedron  is  located.  The 
(Na3)^'^  clusters  are  assumed  to  be  an  eq\ulateral  triangle  characterized  by  a 
distance  d2  from  a  Na  atom  to  the  cell  center.  As  the  surrounding  groups 

become  smaller,  the  electron  ground  state  is  more  extended.  If  d2  =  1.38A,  for  the 
d-model,  the  probability  of  finding  the  electron  within  3.4A  firom  the  tetrahedron's 
center  is  only  3%;  the  probability  of  finding  the  electron  around  the  (Na3)^'^ 
clusters  exceeds  90%.  This  is  inconsistent  with  the  esr  data  which  suggest  that 
the  electron  is  localized  within  the  Na4  cluster. 

V.7.  Polarization  effects 

The  existing  experiments  have  been  performed  with  powders  but 
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experiments  with  doped  single  crystals  are  possible.  In  this  case  one  controls  the 
direction  of  the  electric  field  with  respect  to  the  frame  of  the  zeolite.  If  the  sodium 
clusters  inside  the  zeolite  are  ordered  varying  the  laser  polarization  may  be 
informative.  For  this  reason  we  examine  the  spectrum  for  the  model-c  for  2 
different  polarizations  (Fig.  11).  The  energy  levels  in  the  two  spectra  must  be 
identical  and  the  only  difference  is  in  the  intensities.  The  appearance  of  a  new 
peak  does  not  mean  that  the  system  has  a  new  eigenstate  but  that  the  change  in 
polarization  has  changed  the  oscillator  strength  of  the  state. 

V.8.  Aikafi  substitution  effects 

It  is  not  difficult  to  dope  sodalite  with  other  alkali  atoms.  In  anticipation 
of  such  experiments  we  have  examined  the  absorption  spectrum  obtained  by 
adding  li  or  K  atoms  (Fig.  12).  From  top  to  bottom,  the  three  panels  in  Fig.  12  are 
spectra  of  the  K,  Na  and  li  doped  sodalite,  respectively.  In  our  model  these  atoms 
differ  through  the  cutoff  distance  in  Eq.(in.l)  (see  Table  IV).  The  most 
notable  trend  is  a  blue  shift  in  going  fin)m  K  to  Na  to  Li  impurity. 
Acknowledgements:  The  authors  are  grateful  to  the  Science  and  Technology 
Center  QUEST  and  to  the  Office  of  the  Naval  Research,  for  supporting  this  work. 
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FIGURE  CAPTIONS 

Fig.  1:  One  cage  in  the  sodalite  matrix.  The  A1  and  Si  atoms  are  at  vertices 
and  they  alternate.  There  is  on  0  atom  on  each  line.  The  Na  atoms  sit 
inside  of  the  cages.  The  sodalite  is  a  space-filling  periodic  collection  of  such 
cages. 

Fig.  2:  The  absorption  spectrum  of  the  sodium  atom  for  the  potential  given 
by  Eq.(in.l)  using  x  s  40fs  in  Eq.(IL8).  The  low  intensity  absorption  cross 
section  a(a)) ,  in  arbitrary  units,  is  plotted  versus  energy  in  units  of  eV.  The 
intensities  of  the  first  two  x)eaks,  which  go  off  the  scale,  are  911.8  and  53.9  in 
the  arbitrary  units  used  here. 

Fig.  3:  The  absorption  spectrum  of  the  The  same  parameters  as 

for  Fig.  2.  The  intensities  of  the  first  two  peaks,  which  go  off  the  scale,  are 
1876.5  and  928.3  in  the  arbitrary  units  used  here. 

Fig.  4:  The  orientation  of  the  (Na4)'^  tetrahedron  in  the  coordinate  system 
used.  The  polarization  of  the  laser  electric  field  is  specified  with  respect  to 
the  same  coordinate  system. 

Fig.  5:  The  absorption  spectrum  0(0)  of  the  a,b,c,and  d-models  for 
fiamework  charges  (Table  I).  We  used  t  s  30  fs  in  Eq.(IL8). 

Fig.  6:  A  sequence  of  absorption  spectra  with  the  time  resolutions  a)  tsfi  fs, 
b)  ts  20  fs  and  c)  Xs40  fs. 

Fig.  7:  The  time  evolution  of  the  probability  that  the  electron  in  the 
promoted  state  is  located  at  various  distances  fiam  the  center  of  the  cell  (see 
Eq.V.l):  (a)  0  ^  r  ^  3.4  A  (the  solid  line);  (b)  3.4A  ^  r  ^.4  A  (the  short  dashed 
line  ) ;  and  (c)  5.4  A  ^  r  ^12.0  A  (the  long  dashed  line). 
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Fig.  8:  The  absorption  spectra  for  the  d  model  for  various  tetrahedron  sizes 
di.  (a)  di=2.60lA;  (b)  di=2.50lA;  (c)  dj  =  2.401  A ;  (d)  di=2.30lA. 

Fig,  9:  The  spectra  of  the  d-model  sodalite  for  two  different  orientations  of 
the  (Na3)^'^  groups  in  the  sodalite. 

Fig.  10:  Spectra  of  the  d-model  sodalite  for  three  sizes  of  the  (Na3)^'^ 
clusters,  characterized  by  the  distance  d2  from  the  ions  to  the  center  of  the 
cluster  (which  is  an  equilateral  triangle),  (a)  d2=2.50lA;  (b)  d2-2.40lA;  (c) 
d2=2.30lA. 

Fig.  11:  Spectra  of  the  d-model  sodalite  laser  polarizations,  (a)  The  electric 
field  is  in  the  z-direction  (see  Fig.  4).  (b)  The  electric  field  is  in  the  z- 
direction. 

Fig.  12:  The  spectra  of  sodalite  doped  with  (a)  K  ,  (b)  Na  ,  or  (c)  Li.  We 

f 

assume  the  same  structure  for  the  tetrahedron.  The  only  difference  is  the 
electron  alkali  interaction. 
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TABLE.!: 

The  charges  on  the  zeolite  firamewo^  Qj  is  the  charge  and  is  the  cutoff 
distance  used  for  species  i  in  the  interaction  potential  in  Eq.Gn.l). 


Model 

'  a 

b 

c 

d 

e 

Qnh 

+1.0 

+1.0 

+1.0 

+1.0 

+1.0 

Qsi 

0.0 

+1.0 

+1.5 

+1.9 

+2.5 

Qai 

0.0 

+0.8 

+0.85 

+0.9 

+1.7 

Qo 

-0.25 

-0.7 

-0.8375 

-0.95 

-1.3 

1.725 

1.725 

L725 

L725 

1.725 

R‘si 

— 

0.77 

0.68 

0.60 

OJO 

Rai 

— 

L66 

1.55 

1.50 

LOO 

0.20 

0.20 

0.20 

0.20 

0.20 
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TABLE  n: 

The  ionization  potential  (IP),  the  oscillator  strengths  (i^)  and  the  peak  firequendes 
((Dj)  for  our  Na  atom  model,  compared  to  measured  values. 

Na  Model  Na  Measured* 


IP  (eV) 

5.139 

5.139 

h 

0.975 

0.95-1.0 

0.017 

0.014 

^3 

0.005 

0.002 

0)1  (eV) 

2.22 

2.10 

flJa  (eV) 

3.78 

a75 

0)3  (eV) 

4.59 

4.34 

■  (a)  Measured  values  of  the  oscillator  strengths  are  taken  firom  Ref.20  [  H.G.  Kuhn 
•Atomic  Spectra,  Second  Edition,  Longmans,  London,  1969,  pl73]  while  the  peak 
firequendes  are  taken  firom  Ref  21  [AJL  Stiganov  and  N.S.  Sventitskii,  Tables  of 
Spectral  Lines  of  Neutral  and  Ionized  Atoms,  Plenum,  New  York  ,  1968  ] 


TABLE  m: 

ProbabiKty  distributions  of  the  ground  state  [  H  ,  rj )  from  Eq.CV.D]  for  variotis 
framework  charges  (  models  a,b,c,d  and  e  in  Table  Distances  are  given  in 
Angstroms. 


Model 

a 

b 

e 

d 

e 

P(  0.0, 3.4) 

0.97 

0.97 

0.97 

0.96 

4.(-3) 

P(  3.4, 5.4) 

2.7(-2) 

2.6(-2) 

2.8(-2) 

3.7(-2) 

0.12 

P(  5.4, 12.0) 

6.(-4) 

5.(-4) 

7.(-4) 

2.(-3) 

0.87 
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TABLE-IY: 

The  ionization  potentials  (IP)  and  the  cutoff  distances  used  in  Eq.(III.l)  for 
li,  Na  and  K  impurities. 


Element 

li 

Na 

E 

R*i  (A) 

1.61 

L725 

2.21 

Calculated  IP  (eV) 

5.39 

5.14 

4.34 

Meastnred  IP  (eV) 

5.39 

5.14 

4.34 

arb. 


^aV 


arb.  units 


arb.  units 


arb.  units 
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